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A self-contained mapping closure approximation 

for scalar mixing 

By Guo- Wei He f AND Zi-Fan Zhang + 


1. Objective and motivation 

Scalar turbulence exhibits interplays of coherent structures and random fluctuations 
over a broad range of spatial and temporal scales. This feature necessitates a probabilistic 
description of the scalar dynamics, which can be achieved comprehensively by using prob- 
ability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs 
(Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three 
dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), 
the advection and reaction can be treated exactly but the effect of molecular diffusion 
has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion 
can be expressed as conditional dissipation rates or conditional diffusions. The currently 
used models for the conditional dissipation rates and conditional diffusions (Pope 1991) 
have resisted deduction from the fundamental equations and are unable to yield satisfac- 
tory results for the basic test cases of decaying scalars in isotropic turbulence, although 
they have achieved some success in a variety of individual cases. The recently developed 
mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; 
Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates 
and conditional diffusions, and the models obtained can successfully describe the shape 
relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. 
However, the mapping closure approach is not able to provide the rate at which the scalar 
evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach 
is not closed. In this Letter, we will address this problem. 

The evolution rate of scalar is a key quantity in modeling turbulent mixing for both 
conserve and reactive scalars (Cha & Trouillet 2003). It specifies the characteristic time 
scale of scalar evolution. It has been shown that the decay rate of scalar depends on the 
relative length scale ratio of the initial scalar and velocity fields (Warhaft & Lumley 1978; 
Sreenivasan, Tavoularis, Henry & Corrsin 1980; Durbin 1982; Eswaran & Pope 1988; 
Mell, Kosaly & Riley 1991), and recently, the asymptotic decay of scalar turbulence 
has been extensively studied (Eyink & Xin 2000; Chaves, Eyink, Frisch & Vergassola 
2001; Chertkov & Lebedev 2003). Nearly all existing models for scalar mixing, ranging 
from the simple (conditional) moment approaches to the full PDF approaches, require 
information on the time scales. These models are mainly based on the assumption of 
a direct proportionality between the scalar time scales and the turbulence time scales 
(Fox 1995). Moreover, they exclude the effects of chemical reaction on the time scales 
of scalar evolution. The mapping closure approach of time-dependent reference fields 
(Girimaji 1992) can provide the time scale externally, which highlights an attack line to 
this problem. 
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We develop here a mapping closure approximation (MCA) approach for the time scale 
of scalar evolution. The MCA approach is based on multi-point joint PDFs. In the map- 
ping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989), an unknown random 
field is mapped from a known random Gaussian field so that the evolution of the un- 
known random field can be described by the mapping function. The mapping function 
is obtained from the transport equation for the unknown random field and the Gaussian 
closure. Since the mapping function is constructed at the level of one-point PDFs, it is 
not able to provide the information on two-point statistics, such as the time scales (He, 
Rubinstein & Wang 2002). In the MCA approach, the mapping functions are constructed 
at the levels of multi-point joint PDFs. The mapping function based on the two-point 
joint PDFs could provide the correct information on time scales. 


2. Main results 

We consider the simple case of a reactive scalar advected by a stochastic velocity field: 

^ + u • V<p = TVV + Q(ip), (2.1) 

where the velocity field u obeys V • u = 0 and is independent of the scalar field. With- 
out the loss of generality, it may be prescribed as a known homogeneous and isotropic 
Gaussian field. T is a molecular diffusivity, Q(ip) mimics a one-species chemical reaction. 

In the MCA approach, the scalar field is mapped from a known random field by a 
mapping function 


(2.2) 

Here, the known random field 0(x,t) is taken as a Gaussian reference field. Its one-point 
and two-point joint PDFs are defined by 
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Here r is the magnitude of separation vector r. The mapping function is required to rep- 
resent the one-point and two-point joint PDFs of the scalar ip via the following equations 
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In the classic mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989), 
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the mapping function X is only required to represent the one-point PDF of the scalar ip 
via equation (2.6). Differentiating equation (2.6) with respect to t yields 


dfi 

dt 


d_ 

dip 



= 0. 


( 2 . 8 ) 


Meanwhile, the transport equation for the one-point PDF, /i(^i>,£), can be derived by 
the test function method (Gotoh & Kraichnan 1993; Kimura & Kraichnan 1993) as 


^ + ^[/ 1 < rv V + Q ( ^)IV')]=°- (2-9) 

Thus, comparing the two equations (2.8) and (2.9) with the substitution of (2.2), we 
obtain the exact result 


¥ = r<vVk = %t)) + W 


(2.10) 


The conditional moment in (2.10) can be evaluated from the mapping function (2.2) and 
the Gaussianity (2.3) and (2.4) of the reference field 9 (Panchev 1971) 


r(vV| 'p = x(T),t)) = -c f f'(o,t)r 


d 2 X dX 
V- 


( 2 . 11 ) 


[ dr] 2 dp 

where the prime denotes the derivative with respect to separation r. As a result, the 
transport equation for the mapping function (2.2) becomes 
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It is easily shown from the Gaussianity (2.3) and (2.4) that ((V#) 2 ) = —Cp"(0,t), 
where (7 = 2 for a two-dimensional physical space and (7 = 3 for a three-dimensional 
physical space. Equation (2.12) has been obtained in (Chen, Chen & Kraichnan 1989), 
where — Cp (0, t) is represented by the variance ((Vd) 2 ). However, the correlation p(r,t) 
in equation (2.12) still remains to be unknown and has to be input externally. For ex- 
ample, it is set using the results from direct numerical simulation in (Chen, Chen & 
Kraichnan 1989). Therefore, equation (2.12) is unclosed. 

The two-point correlation p(r, t ) cannot be obtained from the one-point PDF cq (p) . 
Rather, it has to be calculated from the two-point joint PDF g 2 (pi,p 2 ,r,t). Hence, we 
propose to invoke the two-point joint PDF (2.4), which is not used in the classic mapping 
closure approach. By differentiating (2.7) with respect to t, we obtain 
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The transport equation for the two-point joint PDF f 2 derived from the test function 
method (Gotoh & Kraichnan 1993; Kimura & Kraichnan 1993) has the form 


-QjT + Vr • [/ 2 ((u 2 -ui) | = 

-^[/ 2 <rvV + 3(^)1^,^)]- 


7^-[/ 2 (rvV 2 + Q{ip 2 )\ipi,ip 2 )]- (2.14) 



280 


Guo-Wei He & Zi-Fan Zhang 
Subtracting (2.13) from (2.14) leads to 

— • [/ 2 UU 2 - Ui) | 1pl,1p2)] = + ~K~j- 1 / 2 ^ 2 ], (2.15) 

32 ot dtp 1 dtp 2 

where 

H k = r(v V I i’k) - r(vVfc I ^>. (2.16) 

The first term in (2.16) is the conditional diffusion on a single given scalar and is 
calculated in equation (2.11). The second term is the conditional diffusion on two given 
scalars at two different locations and, again, can be evaluated using the mapping function 
(2.2) and the Gaussianity (2.3) and (2.4) of the reference field 9 (Panchev 1971) 
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where l = 1 for k = 2 and ^ = 2 for fc = 1. Multiplying (2.15) by -i/q and tp 2 and 
then taking the mean with substitution of (2.2), (2.4), (2.7) and (2.17), we obtain the 
transport equation for p(r, t) as follows 
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Equations (2.12) and (2.18) form a closed system for the mapping function, where 
equation (2.12) describes the evolution of the shape of the mapping function and (2.18) 
specifies the rate at which the mapping function evolves. In equation (2.18), the second 
term on the left-hand side corresponds to advection, the first three terms on the right- 
hand side correspond to diffusion and the last term on the right-hand side corresponds to 
the effect of nonlinear mapping. The last term vanishes if the mapping function is linear. 
We note that p(r, t) is the correlation function of the reference field and is dependent on 
the mapping function. 

The realizability condition of equation (2.18) is \p(r,t)\ < 1- For the pure diffusion 
processes with the initial Gaussian distributions of positive correlations, the diffusion 
term p (r, <) + p (r, t)/r and the damping term —Cp (0, t)p(r, t) decrease the amplitudes 
of the correlation p{r,t), so that the solution of equation (2.18) is realizable. 

The mapping equation (2.12) is closed using the two-point joint PDF constraint (2.7), 
from which the correlation equation (2.18) is derived. Another possibility for the clo- 
sure is to use the constraint of the joint PDF for the scalar and its derivative (Chen, 
Chen & Kraichnan 1989), which leads to an unclosed equation for p(0,f) and its spatial 
derivatives. It points to another direction to go beyond the one-point mapping (2.2) for 
different purpose. 

The performance of the MCA models (2.12) and (2.18) are evaluated against direct nu- 



A self-contained mapping closure approximation 


281 



Figure 1. Pure diffusion for the conserved scalar with a double-delta initial field. The left: 
variance and dissipation vs time. Solid lines for MCA and squares for DNS; The right: probability 
density function. Solid, dashed and dash-dotted lines for MCA at time 0.05,0.15 and 0.5, and 
circles, squares and deltas for DNS at the same time, respectively. 


merical simulation. The scalar equation (2.1) and the MCA modeled system of equations 
(2.12) and (2.18) are numerically solved in a cyclic square of side 27r, using second- 
order Aclams-Bashforth scheme in time and fourth-order central finite-difference scheme 
in space. In all the cases, the non-dimensional molecular diffusivity T = 0.01. Boundary 
conditions are periodic in space, except that the ones in the direction of the reference field 
are obtained by extrapolation. The initial fields for (2.1) are double-delta distributions 
or isotropic Gaussian distributions with their energy spectra E^k) oc At 1 '/ 3 . Thus, the 
initial mapping for (2.12) and the initial correlation for (2.18) can be calculated from 
their definitions. 

The velocity field is a given homogeneous isotropic Gaussian process with spectrum of 
the form E{k ) oc fc -5 / 3 . We will here present the results for the rapidly changed velocity. 
That is, the velocity fields are set to change at each time step ( Kimura & Kraichnan 
1993). In order to isolate the effects of the MCA models on diffusivity and reaction, the 
advection term in equation (2.18) are calculated directly from the DNS without invoking 
any models. The analytical treatment on this term can be found in ( Kimura & Kraichnan 
1993). 

Figure 1 compares the evolutions of the variance (<£ 2 (x,t)) and the PDF fi(ip,t) ob- 
tained by DNS of equation (2.1) with an integration of the MCA models (2.12) and (2.18) 
for the case of pure diffusion: u = 0 and Q = 0. The initial condition is set as a double 
delta distribution. The results illustrate that the MCA models represent the relaxation of 
the double delta PDF to the Gaussian PDF not only in its shape but also at the correct 
rate of evolutions. 

In Fig. 2, the same plots are made for the diffusion-reaction equation: u = 0 and 
Q(ip) = — 20<p|v?|, with the initial Gaussian distribution. The MCA models are also in 
agreement with the DNS results. It shows that the MCA models can represent the effects 
of both diffusion and reactions. 

Further comparisons are made in Fig. 3 for the advection-diffusion equation with the 
same initial condition as used in Fig. 2, where the Peclet number is about 101. Evidently 
the MCA models are in good agreement with the decay rate of the scalar. 
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Figure 2. Reaction-diffusion for the initial Gaussian scalar. The left: variance and dissipation 
vs time. Solid lines for MCA and squares for DNS; The right: probability density function. Solid, 
dashed and dash-dotted lines for MCA at time 0.0, 0.05 and 0.2, and circles, squares and deltas 
for DNS at the same time, respectively. 



Figure 3. Advection-diffusion for the initial Gaussian scalar. The left: variance and dissipation 
vs time. Solid lines for MCA and squares for DNS; The right: probability density function. Solid, 
dashed and dash-dotted lines for MCA at time 0.0, 0.05 and 0.3, and circles, squares and deltas 
for DNS at the same time, respectively. 


In this Letter, the two-point joint PDF is represented by a known two-point joint 
Gaussian PDF and a mapping function. The correlation of the joint Gaussian PDF and 
the mapping function evolve under their dynamics. Therefore, the representation of the 
two-point joint PDF evolves with the dynamics of equation (2.1) via the correlation 
and mapping function. If it happens that p(r,t) = 0, equation (2.7) is simplified to 
f‘ 2 (id i t) = t)f i(tp 2 , t). This is the stochastic ansatz in the BBGKY hierarchy in 

statistical mechanics (Balescu 1975), where the correlations are specified a priori and then 
fixed. The MCA approach does not invoke this kind of a priori assumption. Moreover, 
a general fV-point (N > 2) joint PDF can also be obtained using the known two-point 
joint Gaussian PDF and the mapping function (2.2), since it can be similarly represented 
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by the mapping function (2.2) and the iV-point joint Gaussian PDF. The latter is fully 
determined by its two-point joint PDF. In this sense, the MCA approach is self-contained. 

3. Conclusion and future work 

A self-contained MCA approach is developed for modeling scalar mixing in a stochas- 
tic velocity field. It differs from the classic mapping closure approach in that the MCA 
approach makes use of two-point PDFs to represent the time-evolving correlations of the 
reference fields and thus the scalar fields. Unlike usual treatments in the BBGKY hierar- 
chy (Balescu 1975), where the representations are specified a priori, the representations in 
the MCA approach are allowed to evolve in coordinate with the dynamics of scalar mixing. 
The results obtained using the new approach are in good agreement with the DNS results 
for the three cases of pure diffusion, diffusion-reaction and advection-diffusion. The ap- 
proach is under further development for more complex situations including multi-scalar 
mixing and inhomogeneous scalar fields, using time-evolving Gaussian or non-Gaussian 
reference fields. 
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